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1. Introduction 

The Ki-Ks mass difference is believed to arise from K^-K'^ mixing via second-order weak 
interaction. It is highly suppressed and very sensitive to the V — A structure of the weak vertex 
in the standard model. We proposed a lattice method to compute this mass difference, including 
both short and long distance effects. Our previous preliminary work on 16^ x 16 x 32 domain wall 
fermion (DWF) lattice shows promising results In this work, we extend our calculation to 
include all operators and various kaon masses on the same lattice. The short distance and long 
distance contributions will be discussed separately. Then we will discuss the operator mixing and 
renormahzation. Finally, the resulting mass differences will be given in physical units. 



2. Setup of this calculation 

The first-order, = 1 effective weak Hamiltonian including four flavors can be written as: 
"'^ = % I V,,V*,,{CiQf +C2Qf), (2.1) 

V ^ q,q'=u,c 

where the Vqqi are CKM matrix elements, C,- are Wilson coefficients and we only include the 
current-current operators, which are defined as: 

Qf' = {siq'j)v-A{qjdi)v-A 
Qf' = {siq\)v-A{qjdj)v-A- 

We neglect the penguin operators in the effective Hamiltonian. This is a good approximation since 
these operators are suppressed by a factor T = |VfrfV,*|/|VKrfVj,* | = 0.0016 in a four flavor theory. 
The essential part of this work is to calculate the four-point con^elator: 



G{tf,t2,h,ti) = (0|r [K\tf)Hy,[t2)Hw{ti)K^{ti)] |0) 



n 

where Nk is the normalization factor for the kaon interpolating operator, tf — tk and — ti should 
be sufficiently large to project onto the kaon state. Since we fix tj and tf in the simulation, this 
correlator depends only on the time separation between the two Hamiltonian |?2 — We will 
refer to this quantity as the unintegrated correlator, which receives contributions from all possible 
intermediate states. 

We can integrate the times t\ and t2 in the unintegrated correlator over a time interval \ta,tb] 
and obtain: 

^=^1 I {0\T{¥\tf)Hy,{t2)Hy,{H)^{tt)]\Q) 

mk-e„ y mk-e„ ] ^^■'^> 

+ "-{f\Hw\m){no\Hw\K'')TA, 
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where T = tt, — ta + I is the integration range, {no) is a state degenerate with kaon if such a state 
exists. We call this amplitude the integrated coiTclator. For sufficiently large T, the exponentially 
decreasing terms become negligible. After the subtraction of the exponentially increasing term and 
a possible quadratic term, the remain term is a linear function of integration range T. The coefficient 
of the linear term gives the finite-volume approximation to AMk up to some normalization factors. 
More details about this method can be found in Ref. 

There are four types of contractions contributing to this correlator as shown in Fig. We 
only include type 1 and type 2 diagrams in this simulation. This simulation is performed on a 
Nf = 2 + 1 flavors, 16^ x 32 x 16 lattice with domain wall fermions, the Iwasaki action, a^^ = 1.73 
Gev, m;t = 421 MeV and jjik = 563 MeV. We use 800 configurations, each separated by 10 time 
units. Two Coulomb gauge wall sources are located at f,- = and tf = 27. The two weak operators 
are introduced in the interval 5 <ti,t2< 22. We calculate the four-point function defined in Eq. 2.3 
for all possible choices of ti and t2. The mass difference is obtained from the slope of the integrated 
correlator plot, and the fitting range is T G [9, 18]. 
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Figure 1: Four types of diagrams that contribute to K^-^ mixing. Only type one and type two diagrams 
are included in this work. The black dots represent a matrix insertion for kaons. the dark circles stand for 
the four-fermion operators. 



3. Short distance contribution 

If there is no charm quark in this calculation, the short distance part of AMk will receive power 
divergent contributions from heavy intermediate states. Such contributions are unphysical lattice 
artifacts. To investigate the divergent chai^acter of this short distance contribution in detail, we 
introduce an artificial position-space cutoff radius R. When we perform the double integration, we 
require the space-time separation between the positions of the two operators to be larger than or 
equal to this cutoff radius : 

^{t2-ti)^ + {x2-xiy->R. (3.1) 

We can plot the mass difference AMk as a function of this cutoff radius R. We use 600 configu- 
rations sepai-ated by 10 time units, with m,^ = 412 MeV and mk = 563 MeV. We only include the 
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operator combination Qi • Qi, i.e. both four quark operators are Qi operators. In Fig. Q we show 
the cutoff dependence of mass difference. The blue curve is a naive uncorrelated two parameter fit: 

AMk\R) = -^ + c, (3.2) 
where b and c are constants. The fitting results shows a power divergent short distance contribution. 




1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

Cutoff Radius 



Figure 2: The mass difference for different values of the cutoff radius R. The blue curve is a naive two- 
parameter fit to a l/R^ behavior. 

We introduce a valence chaim quark to control such unphysical short distance effects. In our 
earlier work [p, we found that GIM mechanism substantially reduce the short distance contribu- 
tion. One might expect that the GIM cancellation would leave us a logarithmic divergence of form 
In(nica) on lattice. However, because of the V — A structure of the weak vertices, the GIM cancel- 
lation is complete, leaving only convergent pieces. Thus if the lattice artifacts associated with large 
ntca and the quenched treatment of the charm quark can be neglected, our four flavors calculation 
should capture all important aspects of AM^. 

In Ref. we were unaware of the absence of a ln{mca) term in the lattice calculation with 
valence charm and we performed an explicit Rome-Southampton style, RI/SMOM subtraction to 
remove it. In fact, the subtraction term reported in Ref. [|T]], performed at a scale jx =1 GeV, was 
zero within errors, consistent with the absence of a true short distance contribution to AM^. 

4. Long distance contribution 

In this section we will discuss the long distance contribution to our calculation of AM/f . The 
integrated correlator receives contribution from both short and long distances. Therefore, in this 
section we examine the unintegrated correlators in Eq . [2.3[ , which depends only on the time sepai^a- 
tion T of two effective Hamiltonians. If the separation T is sufficiently large, the contribution from 
the lightest intermediate state will dominate the signal. There is no vacuum state since we don't 
include disconnected diagrams, so the two lightest states are 71° and n-n states. These two states 
have different parity. We can examine them separately by sepai^ating each left-left A5 = 1 operator 
into parity conserving and violating parts. The results presented in this section are for an average 
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over 800 configurations separated by 10 time units, with a valence light quark mass m/ = 421 MeV 
and various strange quark masses. 

We will discuss the parity-odd channel first. For this case, both weak operators are parity 
conserving which implies that all intermediate states have odd parity. In Figj3|, we plot the uninte- 
grated correlator from operator combination Qi-Qi and the resulting effective mass for the choice 
of kaon mass Mk = 0.4848(8). In the plots of the unintegrated correlators we show both original 
results and the results after the subtraction of the tt" contribution. This subtraction is done using 
the {n^\Qi\K^) matrix element from a three point correlator calculation. Since only the tt" term 
should be present for large time separations, we expect that the results after subtraction should 
be consistent with zero for large Th- In the effective mass plots, we calculate the effective mass 
Mx — Mk from the unintegrated correlators, here Mx is the mass of the intermediate state. For this 
parity conserving case, the lightest state is the pion. The "exact" Mj^ — Mk mass obtained from the 
two point correlator calculation is shown in the plots as a blue horizontal Une which agrees well 
with the computed effective mass. 

Next we will examine the case where parity violating operators appear at both vertices, so the 
intermediate states can only be parity even. In Fig. ^, we present the unintegrated correlators for 
the thi^ee different operator combinations evaluated at a kaon mass Mk = 839 MeV. This kaon mass 
is very close to the energy of two pions at rest, so we expect to find a non-zero plateau at large time 
sepai^ation Th- However, our results ai^e extremely noisy at long distance and we are not able to 
identify such a plateau. The lai^ge noise can be explained as follows. Although the signal should 
come from two-pion intermediate states, we will also have noise, whose large time behavior is like 
that of a single pion intermediate state. Then the signal to noise ratio will fall exponentially for 
lai^ge time separation. The situation here is very similar to what is found for disconnected diagrams. 
We also expect that most of the noise comes from type 1 diagrams, shown in Fig. [I], because the 
topology of type 2 diagrams does not allow a single-pion contribution to their noise. This argument 
is confirmed by plotting the results from type 2 contractions only. If we analyze the type 2 diagrams 
alone, and fit the resulting intermediate state masses the results agree with the two-pion mass very 
well, as seen in the lower right panel of Fig. H. 
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Figure 3: A plot of the unintegrated correlator and resulting effective mass for the parity-conserving operator 
combination Qi Qi and a kaon mass k ~ 839 MeV. In the left-hand plot, the red diamonds and blue squares 
show the result before and after subtraction of the tt" term. In the right-hand plot, the red diamonds are 
effective masses obtained from the unintegrated correlator. The blue horizontal line shows the "exact" value 
of Mji — Mk obtained from the two point correlator calculation. 
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Figure 4: The unintegrated correlators for different, parity-violating operator combinations at a kaon mass 
Mk = 839 MeV. We plot both the full results and the results from type 2 diagrams only. The last plot is the 
fitted intermediate state mass (red diamonds) and the sum of that mass and the kaon mass (blue squares) for 
all choices of kaon masses. The results shown in this last plot are obtained from fitting the type 2 diagrams 
alone. 



5. The Kl-Ks mass difference 

In order to calculate the physical Kl — Ks mass difference, we must connect our four-quark 
lattice operators with the physical AS' = 1 effective weak Hamiltonian Hw- Thus, we must deter- 
mine the Wilson coefficients and normalize the lattice operators in the same scheme in which the 
Wilson coefficients are computed. This procedure can be summarized by : 

Hw = ^ L V,.^, £Cp(M)(l+Ar^™),(Z")y.ef''"(M) 

y^an'=u,c (=1,2 

r (5.1) 

y-^q,q'=u,c (=1,2 

All the operator renormalization and mixing are performed at a scale /i =2.15 GeV. The Wil- 
son coefficients Cj^^(/x) are calculated following the formulas in Ref. |^. The matching matrix 
^^Ri^MS |g provided by Christoph Lehner. The lattice operator mixing matrix is obtain from a 
non-perturbative renormalization calculation [Q]. Combining all three ingredients we can obtain 
the coefficients for the bare lattice operators. These results are given in Tab. [I[ 
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Table 1: The Wilson coefficients, the Rl MS matching matrix, the non-perturbative lat — > RI operator 
renormalization matrix and their final product, all at a scale ji =2.15 GeV shown in columns one through 
four respectively. 



CMS 


^2 


Aril = Ar22 


Ari2 = Ar2i 


^11 — -222 


^12 — -221 




'-2 


-0.2976 


1.1391 


-6.562 X 10-2 


7.521 X IQ-^ 


0.5725 


-0.01412 


-0.1693 


0.6119 



We now combine all these ingredients and determine the mass difference AMk in physical 
units. We present results for four kaon masses ranging from 563 MeV to 839 MeV. In the heavy 
kaon cases, especially the 839 MeV kaon, the two-pion intermediate state will be close to being 
degenerate with the kaon. Thus,we should remove the term in Equation 2.4 arising from this 
degenerate n — n state. However, we are not able to identify the two-pion intermediate state within 
eiTors. This implies that this on-shell, two-pion intermediate state contributes only a small part to 
the mass difference in our calculation. Therefore, for all the results presented in this section, we 
neglect possible contamination from a degenerate two-pion intermediate state. 



Table 2: The contribution of the three operator products evaluated here to the mass difference AMk for four 
different choices of the kaon mass. The final column gives the complete long and short distance contribution 
to AMk expressed in physical units. 



Mk (MeV) 






am|2 


AM^(xlO-'2 MeV) 


563 


6.42(15) 


-2.11(16) 


1.56(9) 


5.12(24) 


707 


8.94(23) 


-3.16(27) 


2.26(14) 


6.92(39) 


775 


10.65(29) 


-3.49(35) 


2.67(18) 


8.08(51) 


839 


12.55(37) 


-3.84(46) 


3.11(24) 


9.31(66) 



In summary, we perform a systematic study of the Kl — Ks mass difference calculation using 
the method of lattice QCD. We only include part of the contractions and use a pion mass of 421 
MeV. Our results range from 5.12(24) x lO^'^ MeV to 9.31(66) x lO^^^ for kaon masses 
between 563 MeV and 839 MeV, while the physical result is 3.48 x lO^'^ MeV. 

The author thank very much all my colleagues in the RBC and UKQCD collaborations for 
valuable discussions and suggestions. Especially thanks to Prof. Norman Christ for detailed in- 
structions and discussions. 
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